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Abstract 



A search for the Higgsstrahlung process e + e~ — > hZ is described, where the neutral 
Higgs boson h is assumed to decay into hadronic final states. In order to be sensitive to a 
broad range of models, the search is performed independent of the flavour content of the 
Higgs boson decay. The analysis is based on e + e~ collision data collected by the OPAL 
detector at energies between 192 GeV and 209 GeV. The search does not reveal any 
significant excess over the Standard Model background prediction. Results are combined 
with previous searches at energies around 91 GeV and at 189 GeV. A limit is set on the 
product of the cross-section and the hadronic branching ratio of the Higgs boson, as a 
function of the Higgs boson mass. Assuming the hZ coupling predicted by the Standard 
Model, and a Higgs boson decaying only into hadronic final states, a lower bound of 
104 GeV/c 2 is set on the mass at the 95% confidence level. 
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1 Introduction 



In the Standard Model (SM) and for masses relevant to the LEP energy range, the Higgs 
boson is predicted to be produced principally by the Higgsstrahlung process e + e~ — > hZ 
and to decay dominantly into the bb channel. This is also the case in large domains 
of the Minimal Supersymmetric Standard Model (MSSM) parameter space (the Higgs 
phenomenology is reviewed, e.g. in Ref. pQ). Most of the searches conducted so far at 
LEP, therefore, tag the b flavour to enhance the Higgs boson signal. 

In other scenarios, however, the decay of the Higgs boson into lighter quark flavours 
or into gluon pairs may be important. Such is the case in general models with two 
Higgs field doublets (2HDM) QUO] or other extended models [3]. In order to be sensitive 
to Higgs bosons predicted by such models, the search described here is based only on 
kinematic selections which are insensitive to the hadron flavour present in the final state. 
Such searches have already been reported by OPAL; these were based on data collected 
at energies close to the Z boson resonance [I] and at a centre-of-mass energy (\/s) of 
189 GeV A similar search has also been reported by ALEPH 6J. 

This paper describes a flavour independent search which is based on OPAL data col- 
lected at centre-of-mass energies between 192 and 209 GeV with an integrated luminosity 
of about 420pb _1 . For the results presented, this search is combined with the earlier 
OPAL searches |U|3]. 
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2 Data sets and Monte Carlo simulation 



The OPAL detector is described in Ref. [7j. The events selected for the analysis have 
to satisfy a set of detector status requirements which ensure that all relevant detector 
elements are active. Events are reconstructed from charged particle tracks observed in 
the central tracking detector and energy deposits ("clusters") in the electromagnetic and 
hadron calorimeters. The tracks and clusters are required to pass a set of quality re- 
quirements jHj- In calculating the visible energies and momenta E vis and P v ; s , either for 
individual jets or for the events, corrections are applied to prevent double-counting of 
the energy attributed to the tracks and to the clusters geometrically associated to the 
tracks [Oj. 

The data sets to which the present analysis applies were collected in 1999 at \fs be- 
tween 192 and 202 GeV and in the year 2000 at y/s between 200 and 209 GeV. After 
the detector status requirements the data sample has an integrated luminosity of ap- 
proximately 420pb _1 . The exact amount varies among the different channels (see Table 

HJ. 

A variety of Monte Carlo samples have been generated to estimate the selection effi- 
ciencies for the Higgs boson signal and for the background processes. In order to cover 
the range of energies of the data, the simulations are performed at fixed values of y/s 
between 192 and 210 GeV and for a set of Higgs boson masses. Spline fits are used to 
calculate the signal efficiencies at intermediate values. 

The Higgsstrahlung process is modelled with the HZHA generator [ID] . Samples of 
1000 to 5000 events were produced at fixed masses, between 30 and 120 GeV/c 2 . The 
Higgs boson is required to decay, either according to the SM, or separately to cc, ssor to 
pairs of gluons. 

The simulated background samples typically have more than 30 times the statistics of 
the collected data. The following event generators are used: KK2F ^1] and PYTHIA 
for the process qqft), grc4f [T3|, KORALW [Hj and EXCALIBUR for the four- 
fermion processes, BHWIDE [TH] for e + e~(7), KORALZ [TJj for and t + t~(j), 

and PHOJET jTHj, HERWIG and VERMASEREN |Uj for hadronic and leptonic two- 
photon processes and for e + e~ — > e + e~77. Hadronisation is modelled with JETSET ^2] 
using parameters described in |21j . The cluster fragmentation model in HERWIG is used 
to study the uncertainties due to quark and gluon jet fragmentation. The Monte Carlo 
samples pass through a detailed simulation of the OPAL detector [22] and are subjected 
to the same analysis procedure as applied to the data. 



3 Analysis 

The search described in this paper addresses the Higgsstrahlung process e + e~ — > hZ. 
The neutral Higgs boson h is assumed to decay into quark pairs of arbitrary flavour or 
into gluon pairs. The following hZ final states (search channels) are therefore considered, 
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depending on the decay of the Z boson: the four-jet channel (Z — > qq), the missing energy 
channel (Z — > z/z/) and the electron, muon and tau channels (Z — > e + e~, and t + t~). 

The analysis assumes that the decay width of the Higgs boson is within the range 
10~ 4 < Th < 1 GeV/c 2 . This ensures that the decay of the Higgs boson occurs within 
about 1 mm of the e + e~ interaction point and that the reconstructed Higgs boson mass 
has a width that is dominated by the experimental resolution, between 2 and 5 GeV/c 2 , 
depending on the search channel. The search strategies are similar to those applied 
by OPAL in the search for the SM Higgs boson [23] (see Ref. [21] for the missing energy 
channel) except that the b-tagging requirements are replaced by more elaborate kinematic 
selections. 

In the searches addressing each of the final states, a preselection is applied first which 
strongly reduces the background while maintaining a high signal detection efficiency The 
preselected events are then submitted to a likelihood test, which discriminates between 
the signal and the two most important background sources, 2-fermion and 4-fermion 
processes. Other background processes, in particular 2-photon events, are negligible after 
the first preselection cuts (see Section 3 of [23])- The likelihood function is constructed 
from reference distributions of a number of discriminating variables which are obtained 
from detailed simulations of the signal and background processes. In the four-jet channel, 
these distributions are obtained from a three-dimensional spline fit to the distributions of 
simulated events where the dimensions are yfs, the hypothetic Higgs boson mass (test- 
mass) and the variable itself. 

Finally, a cut is applied on the value of this likelihood function. For each of the 
search channels, the effect of the preselection and likelihood cuts on the data samples, 
the total background and its contributions, and on the signal detection efficiency for two 
test-masses can be followed through Table ^ 

The signal efficiency is evaluated separately for each of the h — ► bb, cc, ss and gg decay 
hypotheses. At a given test-mass, these efficiencies typically vary by about ±5%. This is 
illustrated in Figure ^a) for the search in the four-jet channel. In deriving flavour inde- 
pendent bounds on the hZ coupling, the smallest of these efficiencies is used; it is obtained 
for h— > gg in the four-jet and tau channels and for h — > cc or bb in the missing energy and 
lepton channels. For the decays into light flavours, h — > uu and dd, the efficiencies are 
slightly higher since the jets are better collimated and because weak semileptonic decays 
are absent; this has been verified explicitly using e + e~ — > ZZ events. These minimal 
efficiencies are shown for all but the four-jet search by the curves in Figure ^b). 

3.1 Search in the four-jet channel 

In the four-jet channel the main background arises from the e + e~ — > WW process. Further 
contributions are from e + e~ — > (Z/7)* — > qq and e + e~ — ► ZZ. The analysis described be- 
low is repeated for fixed test-masses, in steps of 250 MeV/c 2 , between 60 and 120 GeV/c 2 . 
The following preselection is applied: 

1. Events must be identified as multihadronic final states [25:. 
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Cut 


Data 


Total 


qq(7) 


4-fermi. 


Efficiency 


in% 






bkg. 






90 GeV/c 2 100 GeV/c 2 




Four-jet Channel 


Luminosity = 424. 3pb 1 




(1) 


39090 


38831.1 


29929.3 


8322.0 


100 


100 


(2) 


13692 


13648.5 


8602.5 


5012.2 


100 


100 


(3) 


4645 


4504.3 


1077.9 


3418.4 


93 


95 


(4) 


4200 


4038.4 


932.7 


3105.7 


92 


94 


(5) 


3695 


3561.3 


603.2 


2958.1 


90 


91 


(6) 


3594 


3447.2 


581.2 


2866.0 


89 


90 


(7) 


2535 


2399.6 


504.2 


1895.4 


81 


83 


(8) 


2081 


1975.3 


477.2 


1498.1 


78 


80 


(9) 


659 


637.4 


155.8 


481.6 


59 


66 


C 


439 


414.0 


103.8 


136.0 


52 


54 




Missing-energy Channel Luminosity = 420. 9pb 


i 


(1) 


9040 


8524.6 


6063.7 


2382.4 


87 


78 


(2) 


2615 


2391.3 


686.0 


1691.2 


80 


73 


(3) 


2462 


2289.9 


665.4 


1614.6 


77 


73 


(4) 


1635 


1598.4 


110.7 


1487.7 


72 


69 


(5) 


650 


605.4 


48.5 


556.8 


70 


67 


(6) 


298 


291.4 


42.3 


249.1 


65 


62 


C 


123 


133.1 


6.3 


126.6 


45 


48 




Electron Channel 


Luminosity = 422. 3pb 1 




(1) 


18042 


18221.3 


12176.4 


6045.0 


92 


97 


(2) 


558 


538.7 


252.8 


286.1 


75 


78 


(3) 


429 


378.6 


171.0 


207.6 


74 


78 


C 


23 


16.6 


0.2 


16.3 


59 


59 






Muon Channel Luminosity 


= 121, Ipl. 1 




(1) 


18008 


18184.6 


8715.5 


9469.0 


88 


92 


(2) 


505 


477.5 


236.5 


241.0 


77 


81 


(3j 


79 


66.1 


32.6 


33.6 


74 


75 


C 


16 


15 


6.6 


8.4 


64.8 


62.4 






Tau Channel Luminosity = 


--- 409.0pb _i 




(1) 


10417 


10082 


5520.1 


4561.8 


83 


78 


(2) 


1652 


1687.6 


187.0 


1500.9 


62 


61 


(3) 


418 


404.5 


99.5 


305.2 


48 


47 


(4) 


358 


343.1 


96.6 


246.3 


47 


47 


C 


3 


8.8 


0.23 


8.57 


27 


21 



Table 1: Number of events selected in the different search channels after consecutive 
cuts. In each case, the final likelihood cut is denoted by C. The number of events found 
in the data is compared to the expectation from simulations. In the four-jet channel the 
numbers up to and including cut (8) are valid for all test-masses; those for cut (9) and the 
final likelihood cut are given for a test-mass of 100 GeV/c 2 . The last two columns show 
the evolution of the selection efficiencies for Higgs bosons of 90 GeV/c 2 and 100 GeV/c 2 
mass decaying exclusively into hadronic final states at 196 and 206 GeV centre-of-mass 
energy, respectively. 
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2. The effective centre-of-mass energy yU (disregarding initial-state photon radiation, 
see Ref. [25J, is required to exceed 80% of the total centre-of-mass energy. 

3. Events are forced into four jets using the Durham algorithm [2E] and are selected if 
the resolution parameter 2/ 34 is larger than 3 • 10~ 3 . 

4. Each of the jets must contain at least two charged particle tracks to suppress events 
with isolated leptons or photons, like e + e~ — > qq£ + £~. 

5. The matrix element MEqcd for the QCD-induced processes e + e~ — > qqqq and 
e + e~ — > qqgg is calculated |2Ij, approximating the parton momenta by the re- 
constructed jet momenta. The matrix element averaged over all possible flavour 
combinations is required to be within the range —3 < In | MEqcd I < — 1- 

6. The x 2 -probability resulting from a four-constraint (4C) kinematic fit which imposes 
energy and momentum conservation is required to be larger than 10~ 6 . 

7. The four-fermion background is reduced by a cut on the matrix element ME4f of 
the process e + e~ — > qqqq, calculated using EXCALIBUR [T3]. In the calculation 
the parton momenta are approximated by the jet momenta resulting from the 4C 
fit and the matrix element is averaged over all flavour combinations. Its value is 
required to be within the range —8.5 < In |ME4f | < —4.9. 

8. The WW — > hadrons hypothesis is tested in a kinematic 6C fit imposing energy and 
momentum conservation and where the invariant masses of the two jet pairs are 
constrained to the W boson mass. To suppress the WW background, the largest of 
the ^-probabilities, P max (WW), for the three possible jet pairings is required to be 
less than 6.3%. 

9. Finally, for each value of the test-mass, a kinematic fit is performed imposing energy 
and momentum conservation and constraining one dijet mass to the test-mass and 
the other to the Z boson mass. In the fit, the reconstructed Z boson mass is allowed 
to vary within its natural width according to a Gaussian distribution 1 . The largest 
of the ^-probabilities P max (Zh) resulting from the six possible jet assignments to 
the Z and the h bosons is required to exceed 10~ 6 . 



The signal likelihood is constructed using the following 6 variables: (1) the maximum prob- 
ability P max (hZ) of the hZ kinematic fit; (2) the Higgsstrahlung matrix element MEhz [2E] 
for the test-mass considered and for the jet combination which yields P max (hZ); the ra- 
tios (3) MEhz/ME 4 f and (4) MEhz /MEqcd; (5) the difference between the maximum 
and minimum energies of the four jets after the 4C kinematic fit; and (6) P max (WW). 
Distributions of these input variables are presented in Figure El while the likelihood dis- 
tributions for two test-masses are shown in Figure El (a) and (b). Events with a likelihood 
larger than 0.1 are accepted. 

The signal efficiency and residual background rates are affected by the following sys- 
tematic uncertainties: (a) Uncertainties in modelling of the momenta, the angular and 
energy resolutions and the energy scale of the reconstructed jets are less than 2% for both 

1 The sensitivity of the search would be slightly lower if a Breit-Wigner distribution were used. 
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the signal efficiency and the background rate. They have been determined by comparing 
calibration data taken at the Z resonance to the Monte Carlo simulation and transfer- 
ring the observed differences to the simulation of the high energy data, (b) Uncertainties 
in modelling the preselection and likelihood variables are less than 3% for the signal and 

4 — 9% for the background, depending on the test-mass. Weights were applied to the simu- 
lated events such that a x 2 < 1 is obtained when comparing the shapes of the distributions 
from the data and the simulation of the background (for each variable separately). The 
difference of the signal efficiency and background of the weighted and unweighted events 
is considered as the systematic error. It has been explicitly checked that a hypothetical 
signal is not hidden by this procedure, (c) Using alternatively JETSET and HERWIG to 
simulate hadron fragmentation yields a difference of 2 — 13% for the background, (d) The 
cross-section of the four-fermion processes, which dominates the residual background, is 
known to within 2% [3U]. (e) Monte Carlo statistics contribute 1 — 5% for the signal 
and less than 3% for the background. Combining all these effects, the total systematic 
uncertainty amounts to less than 6% for the signal efficiency and 5 — 16% for the residual 
background. 

The number of selected events in the four-jet channel with a likelihood value larger 
than 0.5 is shown in Figure|U(a) for test-masses between 60 and 120 GeV/c 2 . The selected 
data samples for mass hypotheses which differ by less than the mass resolution (of about 

5 GeV/c 2 at high likelihood values) are strongly correlated. For a test-mass of 100 GeV/c 2 , 
439 candidates pass the final likelihood cut of 0.1 while 414 ±53 events are expected from 
background processes and 40 events would be expected from Higgsstrahlung if the hZ 
coupling predicted by the SM is assumed and the Higgs boson decays only into hadronic 
final states. The signal to background ratio becomes more favourable for larger likelihood 
values. 

3.2 Search in the missing energy channel 

Signal events in the missing energy channel are characterised by two hadronic jets and 
a missing mass consistent with the Z boson mass. The dominant backgrounds are four- 
fermion processes, in particular from the semi-leptonic decays e + e~ — > WW — > qqi u, 
and the irreducible process e + e _ — > ZZ — > uuqq. Further contributions are from events 
with particles escaping detection along the beam-pipe, for example from Z boson decays 
accompanied by initial-state photons or the untagged two-photon process e + e~ — > e + e~qq. 
The following preselection is applied: 

1. To reject non-hadronic events, at least 7 charged particle tracks are required. At 
least 20% of all tracks must be of good quality (a minimum number of hits are 
required along the track, see Ref. jH]); this is to reject badly measured events, mainly 
two-photon processes and beam-wall interactions. The total transverse momentum 
p t and the visible mass m vis must satisfy 5 xp t + m vis > y/s/2, and the visible energy 
E vis < 0.8-y/i. The energy measured in the forward detector components [7], which 
cover small polar angles, must be < 2 GeV in the forward calorimeters, < 5 GeV 
in the gamma catcher and < 5 GeV in the silicon-tungsten calorimeter jH2]. The 
overall energy observed in the region | cos#| > 0.9, where 9 denotes the polar angle 
with respect to the electron beam, must not exceed 20% of E vis . 
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2. The missing momentum vector has to point to sensitive parts of the detector, 
|cos# m i SS | < 0.95, and the visible momentum must not have a large component 
along the beam axis, |p* is | < y/s/5. 

3. The tracks and clusters in each event are forced into two jets using the Durham 
algorithm. Events with partially contained jets are rejected by the requirement 

| cos#jct| < 0.95 imposed on each jet. 

4. (Z/7)* — > qq events are suppressed by requesting acop > 5° where the acopla- 
narity angle acop is the deviation of the angle between the two jets in the plane 
perpendicular to the beam axis from 180°. 

5. The missing mass, m miss , must be consistent with the Z boson mass: 50 GeV/c 2 < 
m miss < 130 GeV/c 2 . 

6. Identified semi-leptonic WWdecays with energetic, isolated J32| leptons are dis- 
carded. 

The signal likelihood function is constructed from 5 variables: (1) m miss ; (2) |cos^ miss |; 
(3) max| cos ^jetl i.e. the polar angle of the jet closest to the beam axis; (4) the x 2 ~ 
probability P(1C) of a one-constraint (1C) kinematic fit which imposes energy and mo- 
mentum conservation and constrains the missing mass to the Z boson mass; (5) the angle 
between the missing momentum and the jet with the higher energy: cos#j_ m i SS . The 
distributions of these discriminating variables are shown in Figure 03 and the likelihood 
distribution in Figure El (c). Events with a likelihood larger than 0.4 are selected. The 
Higgs boson mass is reconstructed using the momenta provided by the 1C kinematic fit. 

The number of events passing the likelihood selection is 123 (see Table 1) while 133±11 
events are expected from SM background processes. The most important systematic 
uncertainties [21] are from the modelling of the likelihood input variables and from the 
lepton isolation criteria. The signal efficiencies are affected by a total systematic error 
of 2.9%. The Monte Carlo estimates of the signal efficiencies were reduced by 2.5% 
to account for accidental vetoes due to accelerator-related backgrounds in the forward 
detectors. The reduction factor was determined from randomly triggered events. 

3.3 Searches in the electron and muon channels 

The signal events in the muon and electron channels are expected to have two energetic, 
oppositely charged, isolated leptons and two hadronic jets. The dominant backgrounds 
are e + e~ — > (Z/7)* accompanied by initial state radiation and four-fermion processes, 
mainly from WW and ZZ pairs. The preselection is described in the following: 

1. Events without hadronic jets are rejected by requiring at least 6 charged particle 
tracks. The visible energy E v i S must be larger than 0.6 \fs and the component of the 
total momentum along the beam axis must satisfy |p* is | < £^3— 0.5 y/s. This require- 
ment reduces e + e~ — > (Z/j)*j — > qqj and two-photon processes, e + e~ — > e + e~qq, 
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significantly. All remaining events are forced into four jets using the Durham algo- 
rithm allowing isolated leptons to form low-multiplicity jets. Events are considered 
further if the jet resolution parameter t/ 34 is larger than 10~ 4 . 

2. Two oppositely charged electron or muon candidates must be identified, with ener- 
gies larger than 30 (20) GeV for the higher- (lower-) energy candidate. The energy of 
muon candidates is deduced from the momentum measurement in the central track- 
ing chamber; for electron candidates the energy measured in the electromagnetic 
calorimeter is used. The algorithms to identify muons and electrons are described 
in [HIj and [22], respectively. 

3. The remaining part of the event, after the two lepton candidates are removed, is 
reconstructed as a two-jet event using the Durham algorithm. If the lepton candi- 
dates are muons, a 4C kinematic fit requiring energy and momentum conservation 
is performed to improve the energy and mass resolution of the muon pair; the 
X 2 -probability of the fit must exceed 10~ 5 . For both electron and muon candidate 
events, the invariant mass of the lepton pair is required to be larger than 40 GeV/c 2 . 

The signal likelihood is constructed from five variables in the muon channel and nine 
variables in the electron channel. Those in common are: (1) E vis /^] (2) log l0 yu; (3-4) 
the measured transverse momenta of the two lepton candidates ordered by energy and 
calculated with respect to the nearest jet axis, used to discriminate against semileptonic 
charm or bottom decays; (5) the invariant mass of the lepton pair. For each candidate in 
the electron channel, the additional variables are: (6-7) (E/p— l)/er for the two electron 
candidates, where the momentum p is measured in the central tracking detector, the 
energy E is measured using the calorimeter and a denotes the total error in E/p; (8-9) 
the normalised ionisation energy losses in the central tracking chamber gas [21], for the 
two electron candidates. The event is selected if in the electron case the likelihood is 
larger than 0.3 or in the muon case larger than 0.65. Figures [3] (d) and (e) show the 
distribution of the two likelihood functions. The mass recoiling against the lepton pair is 
taken as the reconstructed Higgs boson mass. 

The number of events passing the likelihood selection is 23 in the electron channel and 
16 in the muon channel (see Table 1) while the corresponding background expectations 
are 16.6 ± 5.1 and 15.0 ± 2.9 events. Systematic uncertainties [21] arise mainly from the 
fragmentation process, determined from a comparison of HERWIG and JETSET, and 
from different four-fermion rate predictions, given by grc4f, KORALW and EXCALIBUR. 
The signal efficiencies have total systematic errors of less than 2% . 

3.4 Search in the tau channel 

Signal events are expected to be composed of two hadronic jets from the Higgs boson 
decay, and two tau leptons from the Z decay. For each of the tau leptons, the decays 
into one or three charged particle tracks ( "prongs" ) are considered, possibly accompanied 
by calorimeter clusters from neutral particles. Important sources of background are the 
processes e + e~ — > ZZ**) — > qq£ + £~, e + e _ — > WW — > qqi^v and e + e~ — > qq(7). The 
following preselection is applied: 
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1. Events must be identified as multihadronic final states The visible energy has 
to exceed 0.3y/s. In order to reject events in which particles escape detection close 
to the beam direction, the missing momentum vector is required to point to sensitive 
detector regions: | cos# m i SS | < 0.95. The scalar sum of the transverse momenta of 
all measured particles has to be larger than 45 GeV/c. 

2. Two isolated tau lepton candidates, each with a momentum between 15 GeV/c 
and 60 GeV/c, are required. These are identified with artificial neural networks 
(ANN) as described in [21] . Separate networks are developed for 1-prong and 3- 
prong decays. From the ANN output, the probability that a candidate is a real tau 
lepton is derived. The probabilities Pi 2 of the two tau candidates are combined to 
a two-tau-likelihood: C TT = P\Pij {P\P<i + (1 — Pi)(l — P2)), which must exceed 0.1. 
If several tau pairs exist in the event, the pair with the largest C TT is chosen. 

3. After removing the two tau candidates, the rest of the event is grouped into two jets 
using the Durham algorithm. A kinematic fit (2C) is applied to the momenta of the 
two tau candidates and the two reconstructed jets, imposing energy and momentum 
conservation. The directions of the tau candidates are approximated by the visible 
momenta of their decay products; their energies are free parameters in the fit. The 
X 2 -probability of the fit must be larger than 10~ 5 . 

4. If both tau decays are classified as 1-prong decays, the momentum sum of both 
charged tracks must be less than 80 GeV/c; this is to reduce backgrounds from 
e + e~ — > TIL — > qq/i + /i~ and qqe + e~. 

The signal likelihood is constructed using nine variables: (1) the visible energy; (2) 
I cos^ m i SS |; (3) 7/34 obtained after reconstructing the event, including the tracks and clus- 
ters of the tau candidates, into four jets (the Durham algorithm is used); (4-5) the angles 
between each of the two tau candidates and the nearest jet; (6) the energy of the most 
energetic muon or electron, if any; (7) the x 2 -probability of a 3C kinematic fit, which in 
comparison to the 2C fit, in addition constrains the invariant mass of the two tau can- 
didates to the Z boson mass; (8) the two-tau likelihood C TT ; (9) the impact parameter 
joint probability of the tau candidate tracks calculated as in Ref. [2H]- The resulting 
likelihood distribution is shown in Figure 01(f). Events with a likelihood larger than 0.8 
are accepted. The invariant mass of the two jets, resulting from the 3C-fit, is taken as 
the reconstructed Higgs boson mass. 

Three events pass the likelihood cut (see Table 1) compared to 8.8±1.5 events expected 
from background. The systematic errors are determined as described in |24j . The largest 
uncertainty arises from the purity of the tau lepton selection. The signal efficiencies are 
affected by a total systematic error of 15 — 17%. 



4 Results 

All search channels combined, 604 candidates are selected, while 588 ± 56 are expected 
from background processes (these numbers apply for a test-mass of 100 GeV/c 2 in the 
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four-jet channel). Figure 0] (b) shows the distribution of the reconstructed Higgs boson 
mass for the candidates selected in the missing energy, electron, muon and tau channels, 
and for the corresponding expected backgrounds, added together. 

No significant excess is observed in any of the search channels over the expected 
background from SM processes. In the following, an upper limit is set on the product 
of the cross-section dhz of the Higgsstrahlung process and the hadronic branching ratio 
Br(h — > hadrons) of the Higgs boson. For this purpose, these search results are combined 
with previous OPAL results, obtained at y/s = 91 GeV in the missing energy, electron 
and muon channels jl], and at y/s = 189 GeV in all channels [5]. 

The limits are obtained by using a weighted event counting method (see Section 5 
of [HB])- The systematic errors are incorporated following Ref. [HI]- The four analyses are 
split into several sub-channels denoted by i depending on the value of the discriminating 
variable. Initially the weight, uJi = Sj / (s, + 6$), is assigned to each channel, where Sj 
and bi are the signal and background expectation, respectively. These weights are a good 
choice if the systematic errors are negligible or do not differ among the channels. To 
obtain weights which are optimal for the general case, the weights are varied until the 
the variance of the distribution of J^i^i ( s « + &i) becomes minimal while keeping J^i^iSi 
constant. The variance is given by 



where I denotes the independent error sources which give rise to the systematic errors on 
the signal, Sjof., and on the background, fejof., of the various channels. A unique solution, 
which minimises the variance, is guaranteed requiring Ui > 0. To fulfil this constraint, an 
unambiguously defined set of sub-channels has to be dropped. 

A test-mass dependent 95% confidence level upper bound k 95 is calculated for the 



where crgj? is the predicted SM cross-section for the Higgsstrahlung process. This bound is 
shown in Figured In calculating this limit, the four-jet and tau channels were considered 
only for masses above 60 GeV/c 2 while the other channels contributed from 30 GeV/c 2 
upwards. Between 12 and 30 GeV/c 2 , only the data taken in the vicinity of y/s = 91 GeV 
are used jl]. The region below 12 GeV/c 2 is covered by a decay mode independent Higgs 
boson search conducted by OPAL [3*3] . 

Limits on the cross-section a for arbitrary Br(h — > hadrons) or for arbitrary hZ cou- 
pling strength g hZ can be derived using the expression 



provided that the effective hZ coupling has the SM Lorentz structure. 

Assuming the hZ coupling predicted by the SM, a Higgs boson decaying only into 
hadronic final states (^(wih) = 1) is excluded for masses up to 104 GeV/c 2 . For a Higgs 
boson also having the decay properties predicted by the SM, this limit is at 100 GeV/c 2 . 




quantity 



k(m h ) 



Chz(^h) x Br(h — > hadrons) 



a 95 = k 95 x <r|^ x Br(h -> hadrons) x (W#hz) 2 ; 
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5 Summary 



A search has been performed for a hypothetical neutral scalar Higgs boson which is pro- 
duced in Higgsstrahlung and which decays to hadrons of arbitrary flavour. The search 
is based on data collected by the OPAL experiment in e + e~ collisions at centre-of-mass 
energies between 192 and 209 GeV. The results have been combined with earlier OPAL 
searches conducted at y/s xs 91 GeV and y/s = 189 GeV. No significant excess has 
been observed over the background expected from Standard Model processes. A mass- 
dependent upper bound is set, at the 95% confidence level, on the product of the Hig- 
gsstrahlung cross-section and the hadronic branching ratio of the Higgs boson. For a Higgs 
boson which couples to the Z boson with Standard Model strength and which decays ex- 
clusively into hadronic final states, a flavour independent lower bound of 104 GeV/c 2 is 
obtained on the mass. 
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Figure 1: Selection efficiencies for the Higgsstrahlung process in the different search channels, 
at y/s = 206 GeV. (a) Four-jet channel, flavour-dependence. The full line shows the result from 
a spline fit to the points with the lowest efficiency, (b) All but the four-jet channel. In each 
case, the lowest of the efficiencies over all hadron flavours is plotted. 
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Figure 2: Distributions of discriminating variables which have been used in the construction of 
the signal likelihood in the four-jet channel, the test-mass being fixed at 100 GeV/c 2 . The 
dots with error bars show the data. The light and dark shaded histograms show the expected 
background from four- and two-fermion processes. The dashed histograms show the signal, scaled 
by a factor ten, expected for a Higgs boson of 100 GeV/c 2 mass, with hZ coupling predicted by 
the SM and decaying only into hadronic final states. 
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Figure 3: Distributions of the signal likelihoods for the searches in the (a-b) four-jet channel, 
(c) missing energy, (d) electron, (e) muon and (f) tau channels. In part (a) the test-mass is 
fixed to 80 GeV/c 2 ; in all other parts it is at 100 GeV/c 2 . The points with error bars represent 
the data. The light and dark shaded histograms show the expected background from four- and 
two-fermion processes. The white histograms added on top of the background contributions 
show the signal expected for a Higgs boson of 100 GeV/c mass (80 GeV/c 2 in part (a)), with 
hZ coupling predicted by the SM and decaying only into hadronic final states. In each case, the 
vertical line indicates the final likelihood cut. 
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Figure 4: (a) Number of candidates selected in the four-jet channel as a function of the test- 
mass mi,, together with the predicted backgrounds and the signal from Higgsstrahlung added 
on top of the background. For the purpose of this figure the likelihood cut is raised to 0.5. (b) 
Combined distributions of the reconstructed Higgs boson mass in the missing energy, electron, 
muon and tau channels. For the signal, the Higgs boson mass is fixed at 100 GeV/c 2 . In both 
parts of the figure, the hZ coupling predicted by the SM and 100% hadronic Higgs boson decays 
are assumed. 
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Figure 5: Distributions of the discriminating variables used to calculate the signal likelihood 
in the missing energy channel. The light and dark shaded histograms show the expected back- 
ground from four- and two-fermion processes. The dashed histograms show the signal, scaled 
by a factor ten, expected for a Higgs boson of 100 GeV/c 2 mass, with hZ coupling predicted by 
the SM and decaying only into hadronic final states. 
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Figure 6: The 95% confidence level upper bound on the product k of the Higgsstrahlung cross- 
section and the hadronic decay branching ratio of the Higgs boson, divided by the Higgsstrahlung 
cross-section in the SM. The thick solid line shows the observed limit. The limit expected on 
average, in a large number of simulated experiments, in the absence of a Higgsstrahlung signal 
is indicated by the dashed line while the dark- and light-shaded areas show the 68% and 95% 
probability bands around the average. 
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